FINITE VOLUME APPROXIMATION OF THE EFFECTIVE 
DIFFUSION MATRIX: THE CASE OF INDEPENDENT BOND 
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Abstract. Consider uniformly elliptic random walk on Z d with independent jump rates 
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, acv< iss nearest . neijihl ><>ur I >< mils of t .he lattice. W e shi >\v th.it the infinite volume effective 

0^ ■ diffusion matrix can be almost surely recovered as the limit of finite volume periodized 

effective diffusion matrices. 



1. Introduction 



Oh 

Consider the following model of random walk in random environment: Given c ^ 1, the 
space of environments 

n = {ab)G[l/c,c],bGZ d *}, 
is specified by the jump rates across the nearest neighbour bonds b E 7L d * . We consider the 
product measure fi = on f2 obtained from a probability measure /in supported on 

[1/c, c\. For every fixed realization £ E f2 of the environment we denote X(t,£) the nearest 
fN| ■ neighbour continuous time random walk on 7j d which jumps according to £(6) rates. The 

symbol E| stands for expectation over this process when the walk starts at x E Z rf . This 
is the quenched regime. 

In the annealed regime the initial random environment is sampled from fi. In the 
uniformly elliptic case we consider here the infinite volume annealed homogenization result 
is well known H,p||,PB,|g|, and the effective diffusion matrix D of X(t) = X(t,£) is 



defined by 



d 



; Vv£R d (Dv,v) = ViD^vj = lim ((X(t, £),«)) 2 . (1.1) 

We shall briefly recall the Kipnis-Varadhan theory and the variational formula for (Dv, v) 
in Subsection |2.1| . Meanwhile let us proceed with describing the finite volume approxima- 
tion for the random environment. This is done in a straightforward fashion: Let N E Z + 
and denote by 7Tjv£ the periodized bond configuration 

(TtfOifc) = & (A) , x E 7L d , x E T N = Z d /2NZ d . (1.2) 

Here x E {— N, ... ,N — l} d , x = x + 2Nz , z E and ttn£ is a periodic configuration 
which coincides with £ on the bonds belonging to the torus T/v- We consider the process 
X(t,7T7v£) describing a particle moving in Z rf according to the rates ttnS,- As we shall 
recall in the beginning of Section || the finite volume homogenization result for X(t, 7Tjv£) 
is, actually, a much simpler statement than the infinite volume one. In particular, there 
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is essentially no difference between the quenched and annealed regimes, and V£ G f2 the 
finite volume effective diffusion matrix Dn(£) is well defined and given by 

d 

Vv G R d (D N (£)v, v)=J2 ViD%te) Vj = lim -E^« ((X(t, n N 0,v)) 2 . 

* — • t — >oo t 

Notice that in general the value of Dn = Dtv(£) depends on the realization of the peri- 
odized environment vr/v£. 

Our main result here asserts that the finite volume effective diffusivities are rapidly 
self-averaging and, furthermore, the sequence {Dn(£)} converges to the annealed infinite 
volume diffusivity D. 

Theorem 1.1. There exists a power v = v(c,d) > 0, such that: 



D%(0-EM(0 >N-n < e~ N \ (1.3) 



for all N sufficiently large. Furthermore, 

lim E^D N (C) = D (1.4) 

Of course, as an immediate consequence of the above theorem we obtain that /i-almost 
surely, 

lim D N (0 = D. (1.5) 

N— >oo 



Our proof of Theorem 1.1 strongly relies on specific properties of the environment mea- 
sure \x\ the concentration type result ( |1.3| ) relies on the independence of the jump rates 
£(6), whereas the proof of ( |1.4| ) makes use of the exchangeability of the joint distribution 
of £(6)-s. Apparently, using the dual variational description of the inverse matrices -D^ 1 
and D~ l , the convergence result (|1.5|) could be established in the general context of er- 



godic jump rates [17]. The finite volume estimate ( |1.3| ) should, of course, require more 
stringent assumptions on the mixing properties of the £-field. Also, as we shall explain in 



Subsection 4.5, our proof has an additional advantage of giving a natural interpretation 
of the translation covariant states for a harmonic interface in a random environment in 
terms of the Funaki-Spohn states ||. 

We would like to mention that approximation results of the type ( |1.4| ) or ([O]) go beyond 
the general spectral analysis of the Kipnis-Varadhan approach as developed in JTo|] or |J. 



This issue is briefly addressed in Subsection ^4j. Similar approximation results have been 



recently derived in [12] in a much more delicate case of the self-diffusion coefficient of the 
tagged particle in the exclusion process. 

Our main motivation for this work came from the theory of massless gradient fields on 
Z d . These are specified by the formal Hamiltonian 

d d 

w) = E E u ( v *^)) = E E u (^)) . 

x i—1 x i=l 

where V, is the discrete lattice gradient; r]i(x) = Vi<fi(x) = <fi(x + ej) — <p(x). In the 
uniformly elliptic case, 1/c ^ U"(-) ^ c, the infinite volume gradient states exist in any 
dimension d ^ 1, and, as it has been established in the paper by Funaki and Spohn |j, 
any translation invariant infinite volume gradient state is decomposable into the convex 
combination of the extremal states (the so called Funaki-Spohn states), which can be 
constructed as limits of finite volume periodized measures with appropriately tilted slopes. 
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The very same approximation by the periodized states leads || to a meaningful definition 
of the (slope-dependent) surface tension. The latter is known to be strictly convex 

On the mesoscopic level, the integral of the surface tension happens to be precisely the 
large deviation rate function for the appropriately scaled height field and it has been 
conjectured in Q that the Hessian of the surface tension governs the equilibrium fluctua- 
tions in the corresponding Funaki-Spohn state. In view of the random walk representation 
developed in EL this conjecture would follow from an approximation result along the lines 



of Theorem 1.1, but for a different type of environment. We refer to |7| for further details. 
To the best of our knowledge this is still an open problem. 

The paper is organized as follows: In Section |2] we briefly recall the variational formula 
for the infinite volume effective diffusion matrix and express it in terms of the corrector 
field. The relevant properties of the corrector field are listed in Proposition |2.1| and, other 
way around, we show that for every ti£l' ! the corrector field ip v is uniquely determined 
by these properties. Subsequently, tp v is recovered as the L 2 (fi) limit of the periodized 



corrector fields ipN,v m Proposition 3.1 of Section ||| This already leads to the convergence 



of the averaged finite volume effective diffusivities E^.D/v(£)> Corollary 3A. In Section || 



we show that Z?at(£) concentrates around its expected value M^Dn^). The concentration 
estimates of Subsection 4.3 prove the main claim in Theorem |1.1| . The key ingredient of 



the proof is an a-priori LP estimate of Subsecti on [4.1| , which is used to derive appropriate 
bounds on the Hamming distance in Subsection [4.2| . By Meyers argument such LP estimate 
follows from discrete version of C alder on- Zygmund inequality. Since we could not find 
direct references to the latter, the proof is sketched in the Appendix. 

2. DlFFUSIVITY IN THE INFINITE MEDIUM 

2.1. The variational formula. We will denote ii(x), i = 1, . . . , d, x G the rate £(&) 
at the bond b = (x, x + e^), with et the unit vector in the i-th. direction. When we write 
only we mean £i(0). T y £, denotes the shifted configuration (r y ^)i(x) = £j(x — y), x,y G Z d , 
i = 1, . . . , d. We also use the notation, for any / : O — > R, 

A/(£) = f(r-e t O - f(t) , D*f(0 = -A/(r ei 6 

In the above notation the associated process of the environment £t = Tx(t)£ seen from the 
particle (see 0, @) is a Markov process with generator 



1=1 

By translation invariance fi is a reversible measure for this process and C is self adjoint in 

It is well known ||] that the annealed effective diffusion matrix D in ( |1 . 1[ ) can be 
recovered from the variational formula: 

d 

(Dv,v) = 2 inf VX + A/) 2 ) (2.1) 

feL 2 (p) ^— j' 

2.2. Corrector field. In general, the above variational problem cannot be solved for /. 
However the next proposition shows that if we only look at "gradients" of / the problem 
has a unique solution. 
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Proposition 2.1. For every v 6 M. d we have 

d d 

(Dv, v)=2Y, fa + = 2 E E ^ fa + V4f (2.2) 

i=l i=l 

where tp l v , i = 1, . . . ,d are the unique elements of L 2 (fi) such that 

• £liA*[&fa + ^(6)] =0, 

• Dkipl = , fj,- a.s. , i,k = 1, ... ,d. 

Proof. For any vector u 6 M rf we define the local drift or current along u as 

= E Wi &(°) " ^" e *)) = " E ^ A^i(0) (2.3) 

i=l i=l 

Let us now define the functions 

X* = (-£ + A)-V», A>0, W G]R d , (2.4) 
or, in terms of the process: 

/•oo 

x£(0=E§ / e" A V,(Tx W e)dt. 



Following Kipnis and Varadhan [10], Theorem 1.3 - see also [16] - one can use the spectral 
resolution of the non- negative operator — C on L 2 (fi) to prove 

XE^f^O, A^O. (2.5) 

Moreover one can prove that there exist functions ip l v £ L 2 (p), i = 1, . . . ,d such that 

d 

5> M &(AxS-^) a ->0, A^O, (2.6) 
i=i 

which in turn by the ellipticity of £ implies strong convergence of DiXv to ip l v in L 2 (/j,). In 
particular, we have 

E^ l<Pv{-£)~ l Vv) = hm E M ^„x« = - lim y^ViE^DiXv = - VV^&V^ • 

A — >0 A — *0 

i=i i=i 

Since ( |2 . 1[ ) is equivalent to 

d 

(Dv, v)= 2J2 M&) v ! " 2E„ (^(-£)" V) , (2-7) 
t=l 



we obtain the first identity in ( |2,2[ ). 

To prove the second one we observe that 



Ee^(^) 2 = -E^ e ^- ( 2 - 



1=1 1=1 
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This follows from fl2.5| ) and the convergence of D%Xv t° 4^1 (P-6D since 
d d 
^E^W) 2 = lim VE^(AX^) 2 = lim E M ( X , A (-£) X 

A — *U A — >U \ 

i=l i=l 



= lim {e^V« " AE M ( X ^) 2 } = - Vi^iti ■ 

Next we show that the functions ip v satisfy the three properties in the statement. The 
first property is obvious since ipl is the L 2 (/u)-limit of a gradient and fi is translation 
invariant. To prove the second identity recall that, in view of (|2.5| ) 

lim E M (/(-£)x$) = Mffv) , f € L 2 (p) . 

A — »U 

Therefore for any / £ L 2 (fi) we have 
d 

E^/A* [6 («< + ^«(0)] = -E^CM) + limE M (/(-£) x A ) =0. 

. A— »U 
1=1 

The third identity follows in a similar fashion from D^DiXv = ^i^kXv^ A > 0. 

Finally, to prove uniqueness we adapt the argument of Theorem 2 in jD| and of Propo- 
sition 3 in (2[: In the language of Q the field ^(0 = (V^(0> • • • > V'uCO) i s > by * ne third of 
the properties in ( |2.2j ), an L 2 (/u)-cocycle. A straightforward modification of the proof of 
Proposition 3 in Q reveals that any L 2 (//)-cocycle . . . , is, in fact, an L 2 -limit 

of gradients: There exists a sequence g 7 £ £ 2 (/u), 7 > 0, with 

E^I^-u*) 2 -^, 7^0, fc = l,...,d. (2.9) 

Uniqueness of ip v is an easy consequence. Indeed, fix v and let u k = ip\ v — ^2„, with 
V'MjV^u both satisfying the required conditions. We then have 

d 

^E^ fc n fc Z) fc / = 0, /£L 2 (/i). (2.10) 
fe=l 



Choosing f = g 7 m ( |2.10| ), with g 7 as in (^), we obtain 

d d 

VE M £ fc (u fc ) 2 = limVE^t/A^ = 0, 

* — • 7 — >U ^ — • 

k=l k=l 

which, by the ellipticity of £, implies u k = fi-a.e., k = 1, . . . , d. □ 



Remark 1. Using the linearity in v implicit in the definition (|2,4| ) we may define the 
functions ip 1 ^ = i\) l e . , i, j = 1, . . . , d, so that 

d 

^ = (^) i = ^^ J - (2.H) 
3=1 

The diffusivity D is then given by the matrix 

D = 2E At A(l + *) (2.12) 

with A denoting the diagonal matrix A iJ = The functions ip 13 are often called the 

corrector fields and the stream matrix. 
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Remark 2. The diffusion constant D has an explicit expression if d = 1. Indeed, for any 
d£R, ^{v + ip v {Cj) must be constant fi-a.s. and normalizing we have 

MO = « [r 1 KrT 1 - 1] , £> = ferT 1 • (2.13) 

By ellipticity we see that the corrector field tp v is uniformly bounded in this case. The 
situation is different if d ^ 2. There is no longer an expression for ip v and all we know 
apriori is that if) € L 2 (p) . For d = 2 this has been upgraded in Q to ^ G L 2+t (fi) for 
some e > 0, which (in two dimensions) leads to an almost sure homogenization result. In 
the sequel, c.f. the remark following Theorem 4.1, we shall establish that ip £ L 2+e (/i) in 
any dimension d 2. 

3. DlFFUSIVITY IN THE PERIODIC MEDIUM 

In the case of the periodized environment ( |1.2j ), the process of the environment seen 
from the particle is equivalent to the projected random walk AT(i,7r/v£) on T/v, i.e. an 
irreducible Markov chain with finite state space and unique invariant measure pn defined 
by 

E pN f = -±- /(*), /:ZV-R. (3.14) 

1 7V| x&T N 

Then, introducing the periodic function 

<Pv ( x ) = <Pv(t-xTTn€) , 
where <p v is defined by (|2.3| ), we arrive as in ( |2.7| ) at the expression 

d 

(D N (0v,v) = 2^E PJV ((^K 2 " 2E PJV (^f(-£iv)-Vf) , (3.15) 

i=l 

where the generator Cn describes the jumps on T/v, i.e. for any periodic function / on Z d 

£-Nf{x) = (vTAr^)i(x) [f(x + Bi) - f{x)\ + (7T N 0i( x ~ e i) [f( x ~ e i) ~ f( x )] • 

Alternatively, identifying periodic functions with functions on T/v we may write 

d 

£nI{x) = - V* ((TJve)i(ar)Vi/(a:)) , / : T^ -» R , 
i=i 

where V are the discrete gradients on the torus T/v • We then have the variational principle 

d 

(v,D N (0v) = 2 inf m J2® PN UmOi(vi + Vif) 2 ) ■ (3.16) 

1=1 

In this finite-dimensional setting the above problem may be solved directly without the 
approximation procedure outlined in the proof of Proposition 2.1. Namely, define 

XN,v( x ) = (-£tf)~V* (s) > ^eZ d (3-17) 
and observe that XNv is a bounded periodic function for each N. Boundedness follows 
from the exponential mixing properties of the process generated by £/v, since E pjv c^ = 0. 
As in the previous section we then have 

d d 

(v,D N (0v) = 2Y^v i E PN {ir N Z) i (v i + ^ v ) = 2 £ E p „ (t^)* (*i + ^,«) * ( 3 - 18 ) 
i=l i=l 
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where tp l N v , i = 1, . . . , d are defined by 

i>N,v{ x ) = Xn,v(x + ti) - Xn,v(x) = ViXN,v(x) (3.19) 
Clearly ifj l N v are periodic functions. The dependence on £ becomes explicit in the notation 

3.1. Convergence. From now on ip l N will be regarded as a function of £, by setting 

fe(0 = ^,*(°'0, i = l,...,d- (3.20) 

Proposition 3.1. For each v G V>at,d — > in L 2 {^l), as N — > oo, where ip v is the 
field defined in Proposition 2.1 

Proof. Notice that, as in (|2i 



]TE pjv (7r,v£M^) 2 = -^^(^Oifc- ( 3 - 21 ) 
i=l i=l 

Then by ellipticity and using Schwarz' inequality we obtain 

maxE pjv (Viv J 2 < <?\v\ 2 . (3.22) 

i ' 

Now we observe that since the bond variables £ are exchangeable under //, the distribution 
of 7Tjv£ coincides with that of t x tt^^ for any x G In particular we see that 

M$N,v? = i^-r E Em(V^(°>™£)) 2 = E M E PiV (^(-,e)) 2 (3.23) 



and ( |3.22| ) provides a uniform bound on the L 2 (/i)-norm of tp l j^ v , which implies weak 
convergence along subsequences. Next we show that any weak limit point must satisfy the 
conditions of Proposition 2.1, thus establishing the weak convergence of ipN,v to t/j v . Let 
„ be a weak limit of ip l N v - 'K ll tp l 00 t, = is a consequence of E,^ip l Nv = 0, which in turn 
follows from exchangeability. To check the second condition we show that for any local 
function / G L 2 {^t) we have 

d 

E E„ [A/6 («i + (0) ] = . (3.24) 

i=l 

From the definition of ip]^ v , writing out ( |3.17| ) explicitly we have 
d 

EVJ ((^(x)^(x,6) = fZ(x), x G T/v . 

i=l 

At x = this yields 

E [( r ^)* ( w < + ^(^tjvO) - 6 («< + fefa^))] = • (3-25) 

i=l 

If / G L 2 (Vt) is local, then /(£) = f(ir N ^) and /(r± ej f) = f(T± ej n N g) when TV is large, 
and by the exchangeability we have 

M(d^)/(TJvO( T «iOA«(TetTi>rO= / M<^)/( r -ei*tfO&fe(f)- 
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Multiplying ( |3,25| ) with / and integrating we see that when N is large 



i=l 



0. 



and the claim ( 3.24| ) follows from weak convergence. The last condition in Proposition 2A 
is proven by similar reasoning, using also 

By the orthogonality relations ( |2.8| ) and (3.21) we, using again exchangeability as in 
fT23D , infer: 

i=l i=l 

In view of the uniform ellipticity condition, £j ^ 1/c > 0, this implies that V'yvu converge 
strongly in L 2 (/j,). □ 

An immediate corollary of the above results is the mean convergence of the diffusion 
constants {Z?at(£)}: 

Corollary 3.1. For every v G M. d , 

lim EpfaDrrffiv) = (v,Dv). 

N—*oo 



Proof. Follows from the representation formulas u2.2|) and ( 3.18 ), exchangeability and the 
L 2 (fi) convergence of ip % N v to tp l v . □ 



4. Almost sure convergence of Dn 

In this section we show that the periodized finite volume effective diffusion matrices 
Dn(<0 converge to the infinite volume effective diffusivity D almost surely in the environ- 
ment, that is (|1.5|) holds. LP bounds on the gradients *^xn,v of the solutions to fl3,17|) set 
up stage for an application of Talagrand's concentration estimates, which imply that for 
every v G 



H-&.B lim \(D N (0v,v)-E tl (D N (0v,v)\ = 0. 

A^oo 



(4.1) 



In the sequel we shall use X% v instead of xn,v, just in order to stress the dependence 
of solutions of ( p. 17 ) on a particular realization £ G [1/c,c] Tn = f2jv of the periodized 
environment. 



4.1. LP estimates. The equation ( |3.17j ) is uniformly elliptic in N G N and £ G O, 



N- 



Theorem 4.1. There exists a power p = p{c) > 2 and a constant a = a(c) < oo, such 
that uniformly in v £R d , N G N and £ G Qn, 



x£T N 



C a\v\. 



(4.2) 
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Remark. By the exchangeability of the ^-environment, ( |4.2[ ) implies 

4>N,v 



supE /t 

N 



P 

< CO, 



where ipN,v has been defined in fl3.20p . Since the infinite volume corrector field could 
be obtained as an almost sure limit of ipN k ,v along a subsequence {N^}, it follows that 
if) v G L p (/j,) for every vel''. 



Because of the uniform ellipticity of jCn the inequality (4.2) follows by Meyers argument 
fll5| (see also Step 1 in Subsection 3.3 in [Q] for the discrete case) once the corresponding 
L q estimate holds for the Poisson equation 

d d 
^V*V,u(x) = -J2V*fi(x) (4.3) 
i=i i=i 

for some q > 2. Namely, we have to show that there exists (5 = /3(g) < co, such that for 
any N G N and every vector field / = ■ ■ ■ , /d) on T\r solutions u of (|4.3[) satisfy: 



11^** ^/%)ll/ft 1<r (4.4) 

( |4.4| ) is a discrete version of Calderon-Zygmund inequality. We could not find a direct 
reference which would cover the case we consider here. For the convenience of the reader, 
a brief sketch of the proof is given in the appendix. 



4.2. Bounds on the Hamming distance. For G Qjy let us define the Hamming 
distance 



d 

dHam.TV = ^ £ ^(aO^G*)' 

xeTjv i=l 



Set q = p/(p — 2) < co, where p = p(c) > 2 is as in the statement on Theorem 4J.. We 
claim that uniformly in v and in £, £' G f^, 



Indeed, by ( gig ) and ( gTlg ), 



M 2 



1 



l^| 



1/9 



(4.5) 



(4.6) 



|Tjv| x G T JvJ =l 
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for every v G M. d and ( G !!jv. Consequently, in view of the a-priori estimate ( [4.2D , we 
infer: 
1 



-\(D N (0v,v)-(D N (£)v,v)\ 



1 n 1 x&T N i=l 



< C 2 |f| 



V? 



i4jEEI«.w-«wi s 



xGTat »=i 



Obviously, 

E E H X ) ~ < C<?d Ham,7V , 

xGTjv i=l 

uniformly in 6 fijv, and the bound fl4.5| ) follows. 



4.3. Concentration of Dn(£). By 

< (Dtv(0v,w) < c|u| 2 , 

for every ti£l li and for each realization of the environment £ € SIat. Thus, given > 
one can find a value Dn{v) £ [0,c|?;| 2 ], such that /^(Ajv(w)) 2eAr/(c|u| 2 ), where we 
define the set An(v) C Oat as: 

^Jv(«) = {^e^iv : |CDj\r(£)«,t;) < cat}. 

Every £' E fi^r such that 

\(D N (e)v,v) -D N (v)\ >2e N , 

certainly satisfies 

min \(D N (£)v,v) - (D N (£')v,v)\ > e N . 

££A N (v) 

By (^^) we arrive to the following bound on such £' in terms of the Hamming distance: 

dHam,7v{^,^7v(w)} = min d Bam>N £} > 04(^5) iV d . 
We are now in position to use the concentration estimate (5.2) in [19]: 



(4.7) 



fM(A N (v)) e [ >| 4 " J 267V 

It remains to choose for some A > 0: Since the sequence of random 

variables (Dn(£)v,v) is uniformly bounded, we readily infer from (4/7) that 

\E fl (D N (C)v,v)-D N (v)\ < 3e N 

for all large enough values of N, and the assertion ( |1.3| ) of Theorem 1.1 follows with 
v = min{A, (d - A)/2q}. 
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4.4. Spectral measures of the local drift <p v . We would like to stipulate the implica- 
tion of (1.5) for the spectral properties of the local drifts ip v as compared to the information 
which could be extracted from the general Kipnis-Varadhan approach: 
In the notation of Section ^ let us introduce the empirical measures 

For every £ G O fixed or, equivalently, for each 7Tjv£ £ fixed there is an obvious 
correspondence between the spaces L 2 (T/v, pjv) and L 2 (Q]y,fj,]y). By ( |2.7] ) and ( p.!5| ) the 
limiting relation (1.5) could be written as 

H-a.s. lim E w (<£„(-£jv)~V«) = e m (W-£)~Vt>) • ( 4 -8) 

Fix now f £ M. d . Let and, respectively, = v% be the spectral measures of <~p v relative 
to the operator (— £,) on L 2 (£l,fi) and, respectively, relative to the operator (— Cn) on 
L 2 (Q/v, /%). Both — £ and all of — £tv are self-adjoint and bounded on the respective 
spaces. Let K be a common upper bound on the spectral radiuses. In terms of spectral 
measures the limit in ( (4.8[) reads as 



N-*coJ r Jo r 

Kipnis-Varadhan approach is based on the fact that (or, respectively, VN ^f r ^ in 

the case of the periodized environment) is bounded above. We claim that our convergence 

result (O) is equivalent to the fi-a.s. uniform integrability of the family |z^(dr)/r|. 

Such equivalence would follow if, for example, we are able to show that /U-almost surely 

the sequence {Hv} weakly* converges to v. The latter is a consequence of 

Lemma 4.1. 

r-K ,-K 

lim / e" nr 4(dr) = / e" n ^(dr), (4.10) 
N ^°° Jo Jo 

fi-a.s for all n £ N. 

Proof. The claim of the lemma essentially follows from the strong law of large numbers: 
Since the local drift (p v is bounded, 

"V w 4(dr) (4.11) 



~ \2(N - M)] d 

L v " xe{-N+M,...,N-M-l} d 

where X(n, ttnS.) is the wrapping round the torus T\r of the random walk X(n,ir^^) 
moving in the periodic environment 7rjv£ and P| denotes the law of such random walk 
with the starting point X(0,ttn^) = x. Since the jump rates of X(t,£) are, uniformly in 
£ £ fl, bounded above by c, 

maxsupPff max \X(t, £) - x\ > M ) ^c 5 e~ C6M2/n . 

x&T N fen X \0 <: t ^ n ' V J 
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Consequently, 



Ei^(X(n,7r N 0) = E|^(X(n,e)) + 0(e 



-c 6 M 2 /n 



(4.12) 



uniformly in x G {— -ZV + M, . . . ,N — M — l} d . Substituting ( [4.12| ) into ( 4.11 ) and choosing 
M = M(N) = y/N, we arrive to the claim of the lemma. □ 

4.5. Some remarks on a massless Gaussian field with bond disorder. An har- 
monic interface in the quenched random environment ^ G is described by the formal 
Hamiltonian 

d 

W) = 2EE6(*)(Vi#*)) 2 - (4-13) 

X 1=1 

We shall see below how the corrector fields introduced in previous sections can be used to 
characterize an interesting class of Gibbs measures for the interaction (4.13). 

In order to localize the interface in d = 1,2 we may consider the field pinned at the 
origin, i.e. we impose (f)(0) = 0. Since the interaction (4.13) is quadratic a full description 
of the set of infinite-volume Gibbs measures for a given typical realization £ G Q, is 
available (||, chap. 13). Namely, it is well known that the set of extremal elements of Q*, 
denoted extC/^, coincides with the set of Gaussian fields on M. zd with a £— harmonic mean 
vector x '■ ^> d ~ * K satisfying 



x(o) = 0, 



(4.14) 



and covariance G^(x, y) given by the Green function of the random walk in the ^-environment 
killed upon hitting the origin. extC/^ is thus characterized by solutions to (4.14). 

To study the tilted states associated to ([4.13 ) it is convenient to work directly with the 

gradient field rji = Vj<^>. This is no loss of information in view of the condition 4>(0) = 0. 
In particular, with the correspondence 0{ = ViX, Q4.14 ) is equivalent to 

^2V*Mx)ai{x)) = 0, V k a i (x) = V i a k {x), i, k = 1, . . . , d , x G Z d , (4.15) 



and extC/^ is characterized by solutions to fl4.15Q , so that v G extQ^ is a Gaussian measure 
on (R rf ) z with mean a satisfying ( 4.15| ) and covariance 



c§(*,v) = 



G^x + e h y + ej) + G^(x, y) - G^{x + e u y) - G^x, y + ej ) . 



Translation covariance and tilted states. A random Gibbs measure is a measurable map 
v : Q — > {£/£, £ G 0} such that G for every £ G 0. The map v is called translation 
covariant when iA o X = u Tx ^ for every £ G fi, x G where 6 X denotes the action 
of translation group and t x is the environment shift. We define the tilted states for our 
interface in random environment as the set of translation covariant random Gibbs measures 



v' such that lA G extQ^. It is not difficult now to use Proposition 2.1 to give a full 
characterization of tilted states. Indeed, since C\- is translation covariant, we only have 
to characterize the translation covariant maps £ — > with obeying ( 4.1 5| ) . For every 
v G M d we know that 

<4(x) = v + ip v (r- x g) 



13 



has the right properties. On the other hand if we additionally require that 0^(0) is in 



L 2 (n,fi) then the uniqueness statement in Proposition 2.1 allows to conclude that erf, 



v £ M. d are the only maps with these properties. In other words, for every v £ M. d the 
Gaussian measure i4 with covariance and mean <rf is the unique tilted state such that 
E^[z4(ry(0))] = v and E At [f|(|ry(0)| 2 )] < oo. The above measures are the analogue in our 
setting of Funaki-Spohn states for translation invariant massless gradient fields (pi). 

Tilted states on the torus and surface tension. To carry this analogy a little further here 
we mimic the construction of ||] to describe tilted states as the infinite volume limit of 
tilted measure on the torus T/v- 

Define Fn as the set of rj £ (M rf ) Tjv such that r\ = V</> for some 4> £ M Tjv , V being as 
usual the discrete gradient on the torus. Consider the probability measure on Tn defined 
by 

D N,v( df l) = ~^H ex P ( - 2 ^2( 7r NOi{.x){fji(x) +Vi) 2 ) m N (dfj) (4.16) 

N,v xGTjv i 

where totv( dfj) stands for the image of Lebesgue measure under the map <j> — * V</>. Notice 
that by symmetry E~« = 0, x £ T\r. Thus when v = 0, v% n is the centered 

Gaussian measure on T^r with covariance - given by the gradients of the Green function 

G^ N (x,y) of the periodized random walk on T\r with killing upon hitting the origin ((4|). 
Now, the linear tilt v only affects the mean in ( f4.16| ) and it remains to compute the mean 
E-5 (fj(x)). We may proceed as follows. 



% W(0))= / (7?(0))dt. 



But 

d - E e (fji(0)) = -cov e (%(0), V V^atO^X^) + 



^ ^2(^NOj( x ) v 3 C N,ij(0,x) = XN,v(ei) ~ XN,v{°) = ^N,v(°i& 

xeT N j 



where in the last line we used the identities (3.17) and ( 3.1S| ) . 



Now for every v £ M. d let i/jL be the probability measure induced by vfr v on the 
field r] = fj + v. The above computation shows that i/fr is the Gaussian measure with 
covariance - and mean vector 

E, (7 7 (x))= W + ^ J v, v (x,0- (4.17) 
From the convergence of C%ij we th en infer, by Proposition 3.1, that at least along a 



subsequence N' — ► oo, v N , ^—almost surely converges weakly to the v— tilted state i/f . 

A last observation is a simple identity connecting the diffusion coefficient to surface 
tension. This type of relation was discussed in p], |?|] in the context of translation invariant 



massless fields. Recalling (4. IE) we define the surface tension on Tjq under tilt v £ K as 



1 , Z N,0 
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By performing the Gaussian integrals and using identities |TL7]) and flgUEp one obtains 
the sought relation 

1 

4 



%(v) = 7 (D N (Ov,v). (4.18) 



An immediate corollary of Theorem LI is then the a.s. convergence of cr^(v) to j(Dv,v). 
As we have already mentioned, the corresponding problem raised in @ in the context of 
an-harmonic models remains open. 

5. Appendix 



The proof of Theorem |4.l| is a straightforward adjustment to the discrete Poisson equa- 
tion (4.3) on Tv of the arguments employed in Chapter 9 of §: Given iV6N and a vector 
field g on Tv define the distribution function rj^g of g as: 

1 

\T N \ 



By the Marcinkiewicz interpolation theorem applied for the solution map T : / h- ► Vu it 
would be enough to check that there exist constants c\ and ci such that for every t ^ 



and, 



W) «S ^^f^, (5.2) 

Since ||Vu||jv,2 ^ 11/11^,2, the second estimate trivially follows from Markov inequality 
(with C2 = 1). It remains to check (|5.1| ) 

To facilitate the exposition we shall consider only the dyadic case N = 2 n . Also, it is 
enough to derive the bound only for the vector fields / of the form f(x) = f(x)e\. 

For every k = 0, 1, . . . ,n — 1 define the decomposition 



Tn = V B * 



x€KT N/K 

of the torus T/v into the lattice boxes B k = x + {0, . . . ,2 k — l} d of the linear size K = 2 k . 
The boxes on different A;-scales are naturally ordered by the inclusion: we say B l y is a 
predecessor of B% if I > k and B k C B l y . By definition T/v is a predecessor of any of the 

boxes B\ on any of the fe-scales. Let us fix a number t > ||/||jv,i- We say that a box B k 
is correct if 

E \m\ = ^ E i/mi < *■ 

By the choice of t the common ancestor T/v is always correct. An incorrect box B k is 
called marked if it is incorrect, but all his predecessors are correct. By construction, 

t < jBj\ E < 2 ^ ( 5 - 3 ) 
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for every marked box B l x . Since all the marked boxes are disjoint, the first of the inequal- 
ities in (|5,3|) implies that the total marked volume 



1 y {Bk{ ^ umi. (5.4) 

\T N \ 1 xl t 

(k,x):B~ is marked 



Let us renumber all the marked boxes as B\, . . . , B[ and decompose / as 

i i 

i=l i=l 
,0 , 



Accordingly, we decompose solutions u of ( [4,3|) as u = u° + Y^i u 3 , where 

d 

^V*V^'(x) = -Vlf(x), (5.5) 

8=1 

for j = 0,1,... ,1. Evidently, 

Since |/°| is bounded above by t, it follows from ( ]5.2| ) that 



. 11/ lllr,2 , ll/°IU,i ,_ „, 

^,V«o(*) < t 2 ~ • (5.6) 

The bulk of the work is, thus, to derive the L 1 estimate on the distribution function 
rj N ^^i ui corresponding to the irregular part Yli f °f the vector field / = fe\. The 

equation ( |5.5D feels the right hand side /* only inside the box Bi. For B, L = B% let us 
define the enlargement Bi = B% via 

B% = = [J By. 

y&KT N / K :\\y-x\\oo <K 

In other words, B\ is the union of B* with all the *-contiguous boxes on the /c-th scale. 
By (0), 

i=l 

In order to estimate r/ N y u , outside Bi write: 

v jU \x) = E ViV^s - y)f(y) = E vh G *( x ~ v)?^) ( 5 - 8 ) 

where 

oo 

and p m is the m-step transition function of the simple random walk on T/v ■ 
There is no loss to assume that f l has zero average: 

E /*(*) = °- 
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Thus, for Bi = By Q , we can rewrite ( fT8| ) as 



V 3 u\x) = Y,(yl j GN{x-y)-Vl J G N {x-y Q ))f\y). 

yeBi 



By Theorem 1.5.5 in |J (d ^ 3) or by Theorem 1.6.5 in |jj (d = 2), 

K 



V hj G N (x - y) - V x jG N (x - y ) ^ c 3 (d)- _ d+1 . 

y& 2/0 1 



It follows that 



for all i = 1, 



E 

xeT N \Bi 

,1. Since the marked boxes B\ 



^ C4 



l/MI, 



E 

x£T N \UBi 



vE 

i=i 



E 



, B[ are disjoint, 



< E i/o»oi- 



As a result, 



1^1 



#{xeT N \UBi : 



vE 



i=i 



> t 



< C 4 



NA 



t 



(5.9) 



which, by (5/7), leads to the desired estimated on the distribution function of the V Yli ^ 
part: 

(4 s II/IIjv,i 



The proof of Theorem 4.1 is concluded. 
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